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Abstract 


Two-phase  continuum  models  have  been  used  to  describe  the  mul¬ 
tiphase  flow  properties  of  solid-gas  and  solid-liquid  mixtures.  The 
approach  is  limited  in  that  it  requires  many  fitting  functions  and  pa¬ 
rameters  to  be  determined  empirically,  and  it  does  not  provide  nat¬ 
ural  explanations  for  some  of  the  qualitative  behavior  of  solid-fluid 
flow.  In  this  report,  we  explore  a  more  recent  single-phase  continuum 
model  proposed  by  Jenkins  and  Savage  (1982)  to  describe  granular 
flow.  Jenkins  and  McTigue  (1989)  have  proposed  a  modified  model 
to  describe  the  flow  of  dense  suspensions,  and,  hence,  many  of  our 
results  can  be  straightforwardly  extended  to  this  flow  regime  as  well. 
The  solid-fluid  mixture  is  treated  as  a  homogeneous,  compressible  fluid 
in  which  the  particle  fluctuations  about  the  mean  flow  are  described  in 
terms  of  an  effective  temperature.  The  particle  collisions  are  treated 
as  inelastic.  After  an  introduction  in  which  we  briefly  comment  on 
the  present  status  of  the  field,  we  describe  the  details  of  the  single¬ 
phase  continuum  model  and  analyze  the  microscopic  and  macroscopic 
flow  conditions  required  for  the  approach  to  be  valid.  We  then  derive 
numerous  qualitative  predictions  which  can  be  empirically  verified  in 
small-scale  experiments:  The  flow  profiles  are  computed  for  simple 
boundary  conditions,  plane  Couette  flow  and  channel  flow.  Segrega¬ 
tion  effects  when  there  are  two  (or  more)  particle  sizes  are  considered. 
The  acoustic  dispersion  relation  is  derived  and  shown  to  predict  that 
granular  flow  is  supersonic.  We  point  out  that  the  analysis  of  flow  in¬ 
stabilities  is  complicated  by  the  finite  compressibility  of  the  solid-fluid 
mixture.  For  example,  the  large  compressibility  leads  to  interchange 
(Rayleigh-Taylor  instabilities)  in  addition  to  the  usual  angular  mo¬ 
mentum  interchange  in  standard  (cylindrical)  Couette  flow. 

We  conclude  by  describing  some  of  the  advantages  and  limitations 
of  experimental  techniques  that  might  be  used  to  test  predictions  for 
solid-fluid  flow. 
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1  INTRODUCTION 


The  flow  properties  of  solid-liquid  and  solid-gas  mixtures  are  important 
for  many  problems  in  science  and  technology,  including  the  properties  of 
molten  lavas  laden  with  rock  crystals,  the  dynamics  of  avalanches  and  rock 
slides,  the  flow  of  red  blood  cells  through  blood,  the  motion  of  sand,  the 
processing  of  paper,  and  the  transportation  of  petroleum  and  coal.  These 
multiphase  mixtures  include  dilute  suspensions,  densely  packed  slurries,  and 
aerated/fluidized  beds.  The  flow  regimes  can  be  homogeneous  or  heteroge¬ 
neous,  laminar  or  turbulent. 

At  the  present  time,  only  rudimentary  progress  has  been  made  in  devel¬ 
oping  a  theory  or  an  experimental  database  for  solid- liquid  and  solid-gas  flow 
even  under  ideal  conditions,  e.g.,  laminar  flow  of  a  monodisperse,  inert,  spher¬ 
ical  particles  through  a  cylindrical  pipe.  Even  if  progress  were  to  be  made, 
the  problems  of  greatest  commercial  interest  —  coal  slurry  transportation 
and  fluidization  of  coal  beds  —  involve  extraordinarily  complex  conditions 
far  from  this  ideal  limit.  For  example,  a  coal  slurry  typically  entails  the  flow 
of  turbulent  water  with  25-50%  concentration  density  of  particles  with  com¬ 
plex,  chemical  interparticle-forces  and  with  a  broad  distribution  of  particle 
sizes  and  shapes.  In  this  sense,  the  realistic  solid-liquid  or  solid-gas  trans¬ 
port  problem  is  much  less  amenable  to  analysis  than,  say,  the  liquid-gas  or 
liquid-liquid  multiphase  flow  problem1  where  both  component  phases  can  be 
simply  described  even  for  realistic  flow  conditions. 

In  developing  an  effective  program  of  research  on  solid  transport,  there¬ 
fore,  a  key  issue  is  the  relative  reliance  on  research  obtained  from  ideal  sys¬ 
tems  versus  efforts  to  improve  empirical  tools  for  realistic  systems.  Clearly, 
empirical  methods  will  remain  the  necessary  and  key  approach  in  engineering 
design,  so  that  the  design  of  new  diagnostic  tools  has  a  high  priority.  For 
example,  the  ultra-sensitive  capacitance  and  fiber-optic  probes  for  measuring 


local,  time-dependent  particle  concentrations  and  velocities,  which  is  being 
developed  by  M.  Louge  and  E.  Giannelis  (Cornell),  appears  to  be  a  technique 
that  might  be  adapted  to  working  conditions  as  well  as  test  beds.  Research 
on  ideal  systems  contributes  insofar  as  it  can  provide  guidance  in  develop¬ 
ing  and  refining  empirical  methods.  For  example,  qualitative  issues  such  as 
the  distinctions  between  flow  regimes,  the  scaling  with  pipe  size,  the  effects 
of  polydispersity,  and  the  role  of  particle-wall  interactions  might  be  usefully 
addressed. 


The  potential  contribution  of  research  on  ideal  systems  depends  largely 
on  the  flow  regime  of  solid  transport,  which  can  be  characterized  by  the 
particle  Reynolds  number,  Re|p: 


Re  Ip  = 


Pol 

Po 


(1-1) 


where  pa  is  the  density,  p0  is  the  viscosity,  and  the  relevant  velocity  is  taken 
to  be  the  local  shear  rate,  7,  times  the  particles  size,  a.  The  cases  in  which 
research  on  ideal  systems  might  best  contribute  are  where  Re\v  is  either  less 
than  one  or  very  high. 


If  Re |p  is  less  than  one,  the  interparticle  dynamics  is  dominated  by  viscous 
hydrodynamic  forces.  The  particle  interactions  can  be  reasonably  treated  by 
simple  kinetic  theory.  Hence,  this  regime  is  the  one  most  reasonably  modelled 
by  ideal  systems,  as  has  been  emphasized  by  Leighton  and  Acrivos2,  and  by 
Brady3  and  coworkers  in  their  studies  of  dense  suspensions.  The  condition, 
Re\p  <  1,  applies,  for  example,  with  slurries  consisting  of  particle  sizes  in 
the  range  a  ~  50/xm  flowing  in  a  six-inch  diameter  pipe  at  5  ft/sec,  in  which 
case  Re |p  «  3  x  10~2.  Even  though  Re |p  may  be  small,  the  Reynolds  number 
based  on  the  macroscopic  flow  can  be  large.  Hence,  one  may  be  able  to  model 
quite  complex  flows. 


For  very  high  particle  Reynolds  number,  we  have  a  different  regime  stud¬ 
ied  first  by  Bagnold.4  In  this  limit,  the  momentum  transfer  and  energy  flow 
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are  both  dominated  by  inelastic  collisions  between  the  particles,  and  the  sus¬ 
pending  fluid  plays  a  secondary  role.  This  regime,  termed  “granular”  flow,  is 
relevant  for  such  geological  processes  as  avalanches  and  sand-dune  formation, 
in  addition  to  its  relevance  for  solid  transport  e.g.,  circulating  fluidized  beds. 
Simple  kinetic  theories  of  particle  interactions  may  be  sufficient  to  explain 
many  phenomena  in  this  regime. 

For  intermediate  ranges  of  particle  Reynolds  number,  progress  is  much 
harder  to  envisage  since  particle  collisions  and  hydrodynamic  forces  compete. 
We  know  of  no  good  theoretical  approach  in  this  regime,  nor  any  experimental 
studies  to  be  recommended.  We  believe  that  attention  should  instead  be 
focussed  on  the  more  tractable  cases  described  above. 

In  addition  to  lying  within  one  of  the  potentially  tractable  regimes,  Re |p  < 

1  or  Re\p  1,  effective  studies  of  ideal  systems  should  conform  to  several 
other  guidelines. 

1.  Particle  properties  should  be  kept  simple.  Effects  such  as  chemical/surface 
interactions,  electrostatic  forces,  van  der  Waals  interactions,  etc.,  should 
be  scrupulously  ignored  until  a  fundamental  understanding  of  simple 
particles  is  obtained.  Fortunately,  experimentalists  appear  to  be  able 
to  design  systems  for  which  this  abstraction  is  close  to  reality.  Acrivos 
and  various  collaborators  have  studied  the  slow  flow  of  suspensions  of 
spheres.1 2 * *  Sayed  and  Savage5  have  focussed  on  granular  flow  under  well- 
controlled  conditions.  In  particular,  they  have  extensively  studied  the 
validity  of  continuum  models  of  granular  flow  and  have  recently  pro¬ 
duced  some  striking  results  concerning  possible  segregation  during  flow 
with  polydisperse  particles  in  both  regimes. 

2.  Hydrodynamic  properties  should  be  tested  in  addition  to  rheological  prop¬ 

erties.  Previous  studies  of  granular  flow  or  the  flow  of  suspensions  have 

emphasized  rheological  properties,  such  as  effective  viscosity  and  shear 


stress.  In  order  to  discriminate  among  theoretical  models,  it  is  im¬ 
portant  to  examine  the  predictions  of  these  models  for  more  complex 
hydrodynamic  phenomena.  For  example,  in  the  theoretical  approach 
we  will  discuss  in  this  paper,  instabilities  such  as  Taylor  vortex  forma¬ 
tion  for  flow  between  rotating  cylinders,  shock  waves  in  grain  flow  past 
blunt  bodies,  and  the  propagation  of  sound  during  shear  flow  could 
be  critical  tests  of  theory.  These  approaches  probably  necessitate  ad¬ 
vances  in  the  ability  to  measure  concentration  profiles  and  correlation 
functions. 

In  this  regard,  we  find  the  new  experimental  efforts  to  measure  particle 
density  and  velocity  profiles  in  solid-liquid  flows  using  refraction-index- 
matched  glass  or  silica-gel  particles  and  transparent  pipelines  to  be  es¬ 
pecially  promising,  this  includes  both  approaches  using  laser  velocime- 
try  (Kadambi,  Case  Western  Reserve)  and  laser  tagging  of  particles 
of  fluid  with  photochemically  sensitive  agents  (Falco,  Michigan  State). 
The  recent  use  of  x-rays  to  study  pattern  formation  in  granular  flow 
(Baxter  and  Behringer,  Duke)  appears  to  be  a  major  breakthrough 
as  well.  These  methods  should  provide  an  important,  new  database 
for  studies  of  monodisperse  and  polydisperse  systems  and  should  be 
continued. 

3.  Large-scale  computation  should  be  limited  to  tests  of  microscopic  prop¬ 
erties  of  multiphase  flow.  Several  excellent  papers  utilizing  molecular 
dynamics  to  study  multiphase  flow  have  appeared  recently.  This  work 
consists  of  directly  modeling  the  interparticle  forces  either  due  to  hy¬ 
drodynamics  (Brady3,  Caltech)  or  to  collisions  (Walton6,  LLNL).  Then 
the  equations  of  motion  can  be  integrated  for  a  small  (50-200)  number 
of  particles  with  boundary  conditions  imposed  to  mimic  a  particular 
flow,  such  as  uniform  shear.  This  work  appears  to  be  well-considered, 
and  is  directly  coupled  to  the  computational  physics  community  so  as 
to  take  advantage  of  advances  in  simulation  capability  (e.g.,  connec- 
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tion  machine  architecture,  cellular  automata,  etc.)  as  they  become 
available.  These  computer  simulations  can  be  used  to  evaluate  the 
rheological  properties  and,  hence,  be  validated  by  the  aforementioned 
experiments.  The  methods  can  also  be  used  to  examine  microscopic 
phenomena,  such  as  the  resuspension  effect  reported  by  Acrivos,  et  al.2 
It  is  clear,  though,  that  this  methodology  will  always  be  limited  to 
small  volumes  and  any  large-scale  flow  is  inherently  untreatable.  For 
that  type  of  problem,  a  less  rigorous  approach,  such  as  a  continuum 
model,  can  be  very  useful. 

4.  For  less  rigorous  approaches,  such  as  continuum  models,  emphasis  should 
be  on  qualitative,  rather  than  quantitative  predictions.  In  some  studies, 
large-scale  computations  and  elaborate  experimental  tests  of  models 
have  proceeded  before  the  basic  assumptions  of  the  model  have  been 
validated.  For  example,  a  leading  theoretical  approach  for  describing 
solid-fluid  flow  is  the  two-phase  continuum  model  (R.  Jackson,  Prince¬ 
ton;  D.  Gidaspow,  IIT).7  This  approach  is  limited  at  the  outset  since  it 
depends  on  many  unknown  functions  that  must  be  fit  to  experiments 
before  any  new  predictions  can  be  made.  More  importantly,  the  sim¬ 
plest  versions  fail  to  properly  predict  some  gross,  qualitative  features, 
such  as  the  sensitivity  of  the  bubble  formation  to  the  viscosity  of  the 
suspending  fluid.  The  response  has  been  to  appeal  to  additional  effects 
requiring  the  fitting  of  new  unknown  functions.  In  this  way,  the  model 
has  evolved  to  be  more  and  more  complicated  with  less  and  less  real 
predictive  power. 

In  spite  of  these  drawbacks,  significant  resources  have  been  focussed 
towards  experimental  measurements  of  the  fitting  functions  and  large- 
scale  computations  of  macroscopic  flow  conditions.  In  one  case,  a 
project  to  develop  high-quality  test  stations  to  study  flow  properties 
in  circulating  fluidized  beds  (T.  Knolton,  IGT)  has  been  slaved,  in 
part,  to  measure  parameters  for  large-scale  computer  simulations  of 
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one  particular  two-phase  continuum  model  (D.  Gidaspow,  IIT).  While 
we  endorse  the  construction  of  the  test  station  and  the  development 
of  an  empirical  database,  we  strongly  recommend  decoupling  the  ex¬ 
perimental  and  theoretical  programs.  In  fact,  we  would  discourage  the 
large-scale  simulations  of  continuum  models  altogether  until  qualitative 
issues  have  been  more  adequately  addressed.  One  needs  first  to  criti¬ 
cally  evaluate  whether  the  entire  approach  is  valid,  perhaps  looking  for 
further  qualitative  tests  of  the  theory. 

As  an  example  of  the  approach  we  are  recommending,  we  present  below 
a  discussion  of  a  promising,  alternative  continuum  model  for  solid-liquid  and 
solid-gas  flow  that  has  been  developed  by  Jenkins  and  Savage,8  McTigue  and 
Jenkins,9  and  HafF10.  The  model  treats  the  solid-fluid  mixture  as  a  single¬ 
phase,  compressible  fluid  in  which  an  effective  temperature  is  introduced 
to  describe  the  particle  fluctuations.  This  model  attracted  our  attention 
because,  compared  to  the  two-phase  models,  it  appears  to  be  a  simpler,  more 
intuitive  approach  with  fewer  unknown  parameters.  We  develop  the  theory 
to  the  point  where  we  can  suggest  simple,  qualitative  tests  to  determine 
whether  the  basic  assumptions  are  workable.  Only  after  such  validation 
would  we  suggest  detailed  computations  of  flow  properties,  such  as  viscosity 
dependence  on  concentration  or  particle  properties.  Also,  we  would  suggest 
a  cautious  apprach  toward  immediate  extensions  of  the  theory  to  the  highly 
turbulent  granular  regime,  as  has  been  suggested  by  M.  Louge  and  J.  Jenkins, 
until  the  model  has  proved  useful  in  the  simpler  flow  regimes. 

The  plan  of  the  report  is  as  follows.  In  Section  2,  we  outline  the  con¬ 
tinuum  theory  of  Jenkins  and  coworkers.  In  Section  3,  we  apply  the  model 
to  various  simple  equilibrium  flow  situations.  We  derive  an  important  con¬ 
sistency  requirement  for  equilibrium  which  amounts  microscopically  to  the 
balance  between  viscous  heating  by  the  shear  and  collisional  dissipation.  We 
also  point  out  some  possible  ways  to  extend  the  current  theory.  For  example, 
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we  introduce  convective-diffusion  terms  to  describe  the  observed  segregation 
in  polydisperse  flows.  In  Section  4,  we  address  the  consequences  of  having  the 
solid-fluid  system  described  by  a  compressible  flow.  We  show  that  granular 
flows  support  sound  waves  which  are  typically  much  slower  than  the  local  flow 
rate.  Hence,  many  systems  will  be  in  the  supersonic  regime  so  that  one  may 
expect  to  observe  shock  waves  of  concentration  in  granular  flow.  (Perhaps 
the  theory  can  be  extended  to  explain  the  concentration  shock  waves  that 
have  been  observed  by  Baxter  and  Behringer  under  conditions  different  from 
simple  shear  flow.)  We  also  discuss  the  consequences  of  fluid  compressibility 
for  the  instabilities,  such  as  in  Couette  flow.  As  an  Appendix,  we  discuss 
existing  and  potential  experimental  imaging  techniques  for  solid-fluid  flow, 
the  critical  component  for  testing  theoretical  predictions. 
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2  SINGLE-PHASE  CONTINUUM  MODEL 
OF  SOLID-FLUID  FLOW 


In  this  section  we  describe  the  single-phase  continuum  model  developed 
by  Jenkins  and  Savage8  to  describe  granular  flow,  a  regime  which  includes 
fluidized  beds  and  rapidly  flowing  sand.  We  will  also  briefly  discuss  more 
recent  work9  by  Jenkins  and  McTigue  using  the  same  framework  to  describe 
dense  suspensions. 


The  model  is  designed  to  describe  a  moderately  concentrated  granular 
fluid  of  monodisperse,  chemically  inert,  spherical  particles.  In  the  next  sec¬ 
tion,  we  will  discuss  how  the  theory  might  be  modified  to  describe  a  fluid 
with  a  distribution  of  particle  sizes.  Moderate  concentration  means  that 
the  particle  diameter,  er,  is  larger  than  the  average  spacing  between  grain 
surfaces,  /t,  but  the  system  is  not  so  densely  packed  that  particles  become 
jammed  together,  thus  preventions  velocity  fluctuations.  To  estimate  the 
relation  between  the  fractional  volume  concentration,  </>,  and  /i/<7,  note  that 
the  fluid  would  become  dense-random-packed  if  each  particle  were  suddenly 
increased  in  diameter  by  h ;  hence, 


(2-1) 


where  <f)onp  ~  -62  is  the  fractional  volume  for  dense-random-packed  spheres. 
By  moderate  concentration,  then,  we  refer  to  4>  between  .2  and  .5  or  he r 
between  .1  and  .5. 


The  model  approach  is  to  treat  the  granular  system  as  a  single-phase, 
compressible,  continuum  fluid.  For  deriving  the  constitutive  relations,  we 
need  to  introduce  the  concept  of  a  “temperature”  related  to  the  local  velocity 
fluctuations  about  the  mean  flow.  If  vT  is  the  average  fluctuating  velocity, 
the  fluctuation  temperature,  T,  is  defined  to  be  ~  Vj.  The  average  time 
between  collisions  is  r  =  — . 

Vt 
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The  model  is  similar  to  the  model  for  a  classical,  dense  gas,  but  there  are 
two  important  distinguishing  features.  First,  we  expect  the  collision  between 
granular  particles  to  be  inelastic.  We  define  a  typical  fractional  energy  lost 
per  collision  to  be  1  —  e;  in  practice  e  ranges  from  nearly  1  (elastic)  to  much 
smaller  values  for  non-spherical  shapes.  In  systems  of  interest,  e  is  probably 
in  the  range  .7-. 9.  Second,  the  source  of  “heat”  in  the  granular  flow  problem 
is  flow  inhomogeneities  or  shear;  there  is  no  direct  analogue  of  heat  sinks  or 
sources  that  can  be  introduced  at  the  boundary. 

The  continuum  equations  follow  straightforwardly  from  the  continuity 
condition  and  conservation  of  momentum  and  energy: 

Tt  =  -'*■*  (2'2> 

dv  -  -*  ~ 

p—  =  -Vp  +  V-t  (2-3) 

2-p^L  =  v  ■  {kVT)  -  pV  ■  v+  7-Vn-r  (2-4) 

where  p  is  the  pressure,  k  is  the  thermal  diffusivity,  and  T  is  the  energy  loss 
rate  per  unit  volume.  The  stress  tensor,  tl3  =  p{di3  +  c^tJV  •  v),  can  be 
expressed  in  terms  of  the  symmetric  strain  tensor,  d,  the  viscosity,  //,  and 
a  constant  of  order  one,  ci.  Note  that  we  intend  the  time-derivatives  to 
be  substantive;  e.g.,  the  left-hand-side  of  Equation  (2-2)  includes  a  spatial 
gradient  of  the  density,  p. 

The  granular  flow  regime  is  defined  as  the  limit  in  which  the  energy 
and  momentum  transport  are  dominated  by  particle  collisions.  Hence,  the 
transport  coefficients,  p,  k,  and  F,  and  the  pressure  can  be  straightforwardly 
estimated  in  terms  of  the  average  fluctuation  velocity,  vT,  the  particle  mass 
m,  the  particle  diameter,  a,  and  the  mean  separation,  h.  It  is  sometimes 
useful  to  express  these  estimates  in  terms  of  the  fluctuation  temperature,  T, 
and  the  particle  mass  density,  p  %  m/<r\ 
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1.  The  viscosity  p  is  proportional  to  the  ratio  of  cross-stream  momentum 
transport  to  the  shear  rate,  7: 


Momentum  flux  = 

Momentum 

Area-Time 

(mAt>)cr-2r_1 

/mcryX  vj 

\  a*  )G  h 

c T 2  . 

~hpVTl' 

Hence,  an  estimate  for  the  viscosity  \s  p  =  p0pT* ,  where  p0  oc  a2 fh. 


2.  Thermal  Diffusivity  is  proportional  to  the  ratio  of  the  energy  flux  to 
the  thermal  gradient: 

Energy  flux  % 

Hence,  the  thermal  diffusivity  is  k  ~  n0pT ?  where  k0  oc  a2/h. 


3.  Energy  loss  rate  per  unit  volume: 


(1  -e) 


mT 

^7 


-  O-Oz-x 


Therefore,  E  =  Y„pT 2  where  T0  oc  (1  —  e)/h. 


4.  Pressure  is  the  force  per  unit  area  produced  by  the  collisions: 

(m&v) 
a2r 

Vf 

PVja  ~h' 

So,  the  equation  of  state  is  similar  to  that  of  an  ideal  gas,  p  =  p0pT , 
where  p0  oc  ajh. 


Pressure 


11 


This  model  is  clearly  simplistic,  even  for  a  monodisperse  distribution  of 
spherical  particles.  Specifically,  it  does  not  take  account  of  rotation  and  slid¬ 
ing  friction  between  particles.  Consequently,  the  model  is  applicable,  at  best, 
over  a  restricted  range  of  concentrations  and  velocities.  Sliding  friction,  for 
instance,  becomes  especially  important  at  low  velocities  and  high  concentra¬ 
tions.  Real  flowing  sand  forms  sandpiles  (or  sand  dunes)  stabilized  by  sliding 
friction  as  the  velocity  becomes  small,  yet  a  continuum  fluid  model  of  the 
type  described  would  predict  a  “sand  puddle”  instead!  Similarly,  Behringer11 
has  noted  some  distinctions  between  flow  of  smooth  versus  irregular  particle 
in  his  studies  of  2D  flow  down  a  hopper.  The  concept  of  local  sliding  is  con¬ 
nected  to  the  local  orientation  of  the  grains  which  has  not  been  included  in 
theories  to  date.  We  can  envisage  a  theory  which  couples  the  translational 
motion  to  the  local  rotation,  which  is  then  hindered  by  friction.  For  now, 
we  will  argue  that  the  present  model  may  be  applicable  for  rather  homoge¬ 
neous  flows  at  intermediate  concentrations  and  high  velocities.  If  the  model 
is  proved  successful  in  such  regimes,  it  will  justify  working  to  enhance  the 
model  to  include  rotation  and  sliding  friction  effects. 


Recently,  Jenkins  and  McTigue9  have  suggested  modifying  this  framework 
to  describe  flows  in  dense  suspensions.  In  this  regime,  particle  fluctuations 
may  play  a  role  similar  to  that  in  granular  flows.  The  key  difference  is  that 
the  interaction  between  forces  is  effected  through  interstitial  viscous  fluid,  via 
lubrication  forces,  rather  than  by  partial  collisions.  Hence  the  dissipation  due 
to  inelastic  collisions  between  particles  (proportional  to  1  —  e)  is  no  longer 
important.  At  high  concentrations,  these  forces  are  dominated  by  lubrication 
forces,  so  F  scales  as  uT/ah ,  where  / 1  is  the  interstitial  fluid  viscosity. 


In  the  viscous  hydrodynamic  regime,  the  moment  um  flux  scales  as  nT’i  /<r, 
(compared  with  pT  for  the  granular  regime).  Hence,  we  must  replace 


P 


(2  -  5) 
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representing  a  change  from  inertial  forces  to  viscous  forces.  By  replacing  p 
in  each  of  the  expressions  above,  we  obtain  the  Jenkins  and  McTigue  model 
for  dense  suspensions. 

One  fine  point  is  that  the  non-Newtonian  character  of  the  flow  (as  ev¬ 
idence  by  a  non-vanishing  contribution  to  the  normal  stress)  is  small  con¬ 
sistent  with  the  fact  that  it  is  strictly  absent  for  spherical  particles  in  the 
hydrodynamic  limit.  In  real  dense  suspensions,  though,  the  normal  stress 
contribution  may  not  be  so  negligible  even  in  the  hydrodynamics  regime,  as 
has  been  stressed  by  Leighton  and  Acrivos.2  This  is  due  to  the  particle  shape 
and  asperity  (surface  roughness)  which  can  effectively  lock  particles  together 
during  a  collision  and  transfer  momentum  perpendicular  to  the  shearing  di¬ 
rection.  It  would  be  useful  to  extend  this  model  (and  the  dynamical  computer 
simulations  of  Brady3)  to  include  this  effect. 


3  LESSONS  FROM  SIMPLE  EQUILIBRIUM 
FLOWS 

In  this  section,  we  apply  the  Jenkins  and  Savage  model  to  various  equi 
librium  flow  problems.  Our  goals  are  to  establish  the  conditions  under  which 
the  model  might  be  valid,  and  to  develop  physical  intuition  about  model 
predictions  for  simple  flow  conditions. 

EXAMPLE  1:  Let  us  suppose  that  we  can  fix  the  entire  mean  velocity 
field  to  be  oriented  along  the  x-direction  with  uniform  gradient  in  the  z- 
direction: 

u  =  uo^l  +  -0x;  (3-1) 

or,  equivalently  there  is  a  constant  shear  rate,  7  =  v0/ L,  where  L  is  the  size 
of  the  system.  In  equilibrium,  the  energy  balance  equation,  Equation  (2-4), 
becomes: 

0  =  V  •  {kXT)+  7  -Vv  -  T  (3-2) 

=  |(^|)+^r-(i)!-r“',rl-  (3'3) 

The  solution  of  this  equation  (also  consistent  with  the  continuity  and  mo¬ 
mentum  balance  equations)  is: 


The  solution  means  that  the  assume  flow  profile  corresponds  to  thermal  equi¬ 
librium  for  fluctuations  in  the  granular  regime.  The  temperature  results  from 
shear:  the  greater  the  shear  rate  7,  the  higher  the  fluctuation  temperature, 
T.  The  single-phase  continuum  approach,  if  it  is  valid  in  any  regime,  should 
be  able  to  predict  properties  of  flows  near  this  limit. 
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Oae  should  also  bear  in  mind  that  any  continuum  model  assumes  that 
the  particle  diameter  a  is  small  compared  to  the  system  size  L.  Hence,  from 
Equation  (3-6),  we  discover  that  the  mean  fluctuation  velocity  vt  is  much 
less  than  mean  flow  velocity  v0.  We  shall  find  this  result  to  be  of  critical 
importance  in  the  sound  speed  analysis  presented  in  the  next  section. 

The  momentum  balance  equation,  Equation  (2-3),  for  the  z  direction 
yields  the  requirement  of  constant  pressure,  which,  combined  with  the  con¬ 
stant  temperature  derived  above,  yields  a  constant  density  profile.  Given  the 
result  for  v r,  this  predicts  that  the  pressure  (normal  stress)  and  the  shear 
stress  ~  /.t are  both  proportional  to  (~))2.  This  prediction  may  be  regarded 
as  a  critical  qualitative  test  of  whether  a  given  granular  flow  is  in  the  range 
of  velocity  and  concentration  parameters  appropriate  for  description  by  the 
single-phase  continuum  model. 

As  it  turns  out,  this  proportionality  to  the  square  of  shear  rate  agrees  over 
the  range  of  moderate  concentrations  and  high  velocities  (perhaps  appropri¬ 
ate  for  fluidized  beds)  in  experiments  by  BagnohE  and  Saved  and  Savage5, 
as  well  as  computer  simulations  of  Walton  et  al..6  (where  sliding  friction  is 
neglected).  Hence,  even  though  we  argued  previously  that  the  model  fails  at 
low  velocities  and  very  high  concentrations  (e.g.,  sandpiles),  there  seems  to 
be  a  physically  interesting  regime  of  concentrations  and  velocities  where  the 
single-phase  continuum  model  may  be  valid. 

EXAMPLE  2:  We  now  give  a  more  precise  treatment  of  the  shear  flow 
case.  We  consider  flow  between  two  parallel  plates  at  ~  =  ±//  at  which 
there  are  imposed  velocities  u(^  =  ±//)  =  ±U .  We  will  assume  that  the 
granular  fluid  obeys  a  non-slip  boundary  condition  -  this  is  almost  certainly 
not  precisely  true  but  our  interest  in  qualitative  findings  allows  us  this  sim¬ 
plification. 

Unlike  example  1  where  we  imagined  fixing  to  be  constant  throughout 


16 


the  bulk,  here  v  is  fixed  only  on  the  boundary;  the  interior  behavior  must 
be  obtained  by  solving  the  equations  of  motion.  By  symmetry,  the  physical 
variables  v,T,  and  p  should  depend  on  z  only,  and  this  dependence  should  be 
determined  by  solving  the  continuity,  momentum  balance  and  energy  balance 
equations,  Equations  (2-2)  through  (2-4).  The  continuity  equation,  Equa¬ 
tion  (2-2),  is  satisfied  trivially  since  v  is  divergenceless.  There  remains  the 
momentum  balance  (along  two  independent  directions)  and  energy  balance 
equations: 

°  =  1(^1)  (3-?> 

0  =  —(pT)  (3-8) 

°  =  "4(',Tif)-^Tf+Mr5(l)2'  (3-9) 

These  equations  imply  T  —  y  with  pQ  a  constant  of  integration  which  must 
♦  be  fixed  by  the  total  density,  and  ~  =  AT ^  where  A  is  a  second  constant  of 
integration. 


Substituting  into  the  energy  balance  equation,  we  find 


d  (  > 

k°t*  K 1 


dT' 
dz  . 


a 2  j 

=  2  ^(KoPoT^)  =  A  K0p0T?, 


where 


(3-10) 


(3-11) 


This  equation  is  directty  solvable  for  the  temperature  field  and  implies 
that  the  temperature  varies,  with  spatial  rate  ^/|Aj.  But  A  is  not  yet  known, 
since  it  depends  explicitly  on  the  constant  of  integration  A.  To  fix  A,  we  need 
to  specify  a  boundary  condition  on  dT jdz  corresponding  to  the  energy  flux 
from  the  walls.  If  dT / dz  =  0,  we  then  must  have 


/l2  = 


r. 

»0 
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corresponding  to  precise  balance  between  energy  loss  and  viscous  “heating” 
due  to  the  shear  flow.  More  generally,  however,  the  flux  will  not  be  zero  and 
in  typical  computer  experiments  appears  to  be  positive  (see  Campbell^12)). 
This  means  that  a  solution  to  the  temperature  equation  has  the  form 

-  E  cosh  Ac  (3-12) 


with  E  determined  via  the  equation 

=  ‘IE2  cosh  MI  sinhA // 


or 

dz 


z-U 


once  A  is  known.  We  then  use  =  AT to  find  the  velocity  field 


v(z)  = 


1/tsinh  \z 


Imposing  v{±FI)  =  ±U  gives  the  relationship 

AE 

(  =  — —  sinh  \Il 


(3-13) 


(3-M) 


determining  A  in  terms  of  A.  This  is  then  combined  with  Equation  (3-11)  to 
find  an  equation  for  A,  so  the  system  closes  consistently. 


The  preceding  analysis  is  important  since  there  appears  to  be  some  con¬ 
fusion  in  the  literature  regarding  exactly  how  many  boundary  conditions 
one  should  impose.  Our  approach  seems  to  agree?  with  that  of  Johnson  and 
Jackson13  but  disagrees  with  that  of  Hanes,  Jenkins  and  Richman14.  In  the 
latter  paper,  there  is  an  additional  restriction  placed  on  the  normal  stress, 
as  opposed  to  considering  just  energy  and  tangential  stress  (which  we  have 
approximated  as  no  slip).  This  leads  to  a  condition  on  the  depth  of  the  layer 
//  which  is  not  physical.  This  is  recognized  by  the  authors  who  relax  the 
density  equation  near  the  boundary.  Our  perspective  is  that  if  the  idea  of 
bulk  flow  equation  makes  sense,  one  must  be  able  to  coarse  grain  over  all 
boundary  effects  and  recover  two  boundary  conditions. 


In  the  next  section,  we  will  discuss  approximate  methods  for  including 
friction  in  the  above  calculation.  For  the  moment,  we  note  that  it  is  clear 


that  for  any  finite  A,  the  above  theory  must  break  down  if  the  depth  H  is 
made  too  large.  This  is  simply  because  the  temperature  decreases  rapidly  as 
we  leave  the  boundary  and  the  density  increases;  at  some  point  the  density 
is  so  large  as  to  cause  the  grains  to  interlock  and  prevent  motion.  In  fact, 
one  item  which  a  full  theory  should  predict  is  exactly  this  maximum  depth 
H.  We  will  return  to  this  in  the  next  section. 


EXAMPLE  3:  Example  2  can  be  modified  to  describe  channel  flow  be¬ 
tween  two  stationary  parallel  plates.  We  choose  the  boundary  conditions 
v(z  =  ±H)  ~  0,  assuming  that  the  velocity  is  nearly  zero  at  the  particle- 
plate  interface.  A  constant  force  F  is  then  imposed  along  the  flow  direction 
(x-axis).  The  equation  so  motion  are  the  same  as  in  Example  2,  except  that 
F  is  added  to  the  right-hand-side  of  the  momentum  balance  equation  for  the 
x-direction.  Consequently,  we  still  have  T  =  ^  (where  pQ  is  an  integration 
constant),  but  now 


p°pT^^A-Fz. 

Because  the  channel  flow  problem  is  symmetrical  about  2  =  0,  there  is  the 


constraint  that 


j=0 


0,  or  A  =  0.  Hence, 


dv 

d~z 


,PoPo , 


T*z. 


(3  -  15) 


If  T  were  constant,  the  solution  would  describe  the  parabolic  velocity  profile 
characteristic  of  channel  flow.  Substituting  this  expression  for  the  velocity 
gradient  into  the  energy  balance  equation,  Equation  (2-4),  we  obtain  the 
same  equation  for  the  temperature  profile  as  before,  Equation  (3-10),  except 
that  A  is  replaced  by: 


A  = 


Vo 

Kq 


KoPoPo, 


(3-16) 


To  describe  the  solution  to  the  temperature  profile  equation  it  is  useful 
to  recast  the  equation  in  the  form: 

a Pe 


dz 2 


+ 


a; 


/Jj  e  =  o, 


(3-17) 
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where  ©  =  T?,  a  =  (—~2  and  8  =  The  general  solution  can  be 

\*oMoP  /  ^ 

expressed  in  terms  of  parabolic  cylinder  functions.1  If  Dv{x)  is  the  parabolic 
cylinder  function  for  real  v ,  then  the  general  solution  to  Equation  (3-15)  can 
be  expressed  in  terms  of  the  real  and  imaginary  parts  of  D~r(etir^*(4a)1^4z), 
where  7  =  i/3(ia)~1^4 .  Linear  combinations  of  these  functions  should  be 
taken  which  respect  the  symmetry  under  2  — *  —z. 

Just  as  in  the  previous  case,  one  needs  to  know  the  boundary  condition  to 
find  the  actual  profile.  Here,  fixing  dT jdz  will,  together  with  the  symmetry 
of  the  temperature  profile,  enable  a  complete  determination  of  0;  unlike  the 
shear  flow,  the  temperature  cannot  be  constant  since  A  is  fixed  by  the  external 
forcing,  not  by  the  internal  dynamics.  All  in  all,  the  lesson  is  that  granular 
channel  flow  is  complicated  compared  to  the  standard  Poiseulle  problem. 

EXAMPLE  4:  In  this  example,  we  return  to  the  problem  of  shear  flow  of 
Example  2,  but  here  with  a  gravitational  field  superimposed  in  the  z  direction 
(perpendicular  to  the  shear  planes).  The  momentum  balance  equation  for 
the  ^-direction  is  changed  to: 

r\ 

g:(pT)  =  -pG.  (3  -  18) 

so  pressure,  p  oc  p7\  decreases  upward.  The  equation  can  be  satisfied  by 
either  a  temperature  gradient  (colder  upward)  or  a  density  gradient  (less 
dense  upward)  or  both.  However,  if  the  continuum  model  is  to  be  valid, 
gravity  must  be  a  small  perturbation  compared  with  the  fluctuation  energy: 
gz  <<  T .  Erom  the  previous  example,  we  know  that  f>  and  T  vary  nearly 
inversely  as  g  — *  0.  Since  the  temperature  is  more  effective  at  balancing  the 
gravitational  force  in  Equation  (3-18)  (the  density  variation  adds  a  term  of 
the  form  ),  the  solution  must  be  that  the  temperature  gradient  will 

point  downward. 
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4  EXTENSIONS 


There  are  several  directions  in  which  it  might  be  possible  to  extend  the 
considerations  of  the  previous  sections.  Perhaps  the  most  crucial  issue  is  how 
to  include  frictional  forces,  but  experiments  involving  particle  size  segregation 
also  may  be  worth  examining.  In  this  section  we  take  a  preliminary  look  at 
these  two  issues.  In  the  next  section,  we  discuss  a  different  direction,  that  of 
acoustic  waves  in  granular  flows. 


Let  us  start  with  the  issue  of  polydispersity.  In  many  experiments  (see  e.g. 
Savage  and  Lun15),  there  is  an  obvious  tendency  of  particles  to  segregate  with 
the  small  ones  falling  to  the  bottom.  Microscopically  this  happens  because 
the  small  particles  can  fall  through  the  holes  made  by  the  dynamical  motion 
of  the  large  particles  but  not  vice  versa.  We  will  try  to  study  a  simpler  case, 
that  of  a  few  large  particle  moving  in  a  continuum  of  small  ones.  We  will 
postulate  that  the  large  particle  concentration  obeys  a  diffusion  equation 

V-  (£>e//Vc)  ~  (4-1) 

We  might  suppose  that  the  effective  diffusivity,  in  analogy  with  Einstein 
formula,  should  be  of  order 


Dcff 


pkT 

6irp.Ro 


for  particles  of  radius  R0  in  a  medium  with  effective  viscosity  p. 


(4-2) 


Recall,  then,  the  case  of  shear  flow  in  a  gravitational  field,  Example  4.  The 
viscosity  is  pT? ,  where  the  temperature  gradient  points  downward.  Hence, 
Dejf  oc  Ti  and  VDejf  is  in  the  —  z  direction.  From  the  above  equation, 

—  —  VDejf  •  Vc  =  De//V2c  (4  —  3) 

implies  that  acts  as  a  negative  convective  term.  The  conclusion  is 

that  the  large  particles  tend  to  be  carried  upward.  In  a  centrifugal  flow,  this 
would  be  upward  and  inward. 
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This  may  be  a  possible  macroscopic  explanation  for  vertical  segregation 
effects.5  Note  that  this  approach  is  consistent  with  the  aforementioned  micro¬ 
scopic  explanation:  all  particles  try  to  move  in  the  forcing  direction  (down¬ 
ward)  but  the  smaller  particles  moves  more  efficiently,  creating  a  new  back- 
flow  of  larger  particles.  Of  course,  since  we  have  treated  the  small  grains 
as  a  continuum  fluid  which  suspends  the  larger  particles,  the  approximations 
breakdown  if  the  larger  (impurity)  particles  were  to  be  replaced  by  impurities 
smaller  than  the  grains.  This  break  down  is  just  as  well  since,  in  this  case, 
the  impurities  should  fall,  not  rise. 

As  the  size  and  density  of  the  impurities  becomes  comparable  to  smaller 
ones,  one  should  take  account  of  their  velocity  fluctuations  as  well.  If  the 
larger  particles  were  to  begin  near  the  bottom  of  the  shear  flow  above,  they 
would  be  in  a  high  T  region,  driven  by  both  a  thermal  and  concentration 
diffusivity  whose  values  could  be  quite  different.  The  situation  is  analogous  to 
double  diffusive  convection  of  salt  in  seawater  heated  by  the  sun.  (Of  course, 
the  geometry  is  upside-down  compared  to  the  granular  flow  problem  in  the 
sense  that  the  higher  temperature  and  concentration  is  at  the  top  for  the 
salt  example  but  at  ihe  bottom  for  the  shear  flow  problem.)  By  this  analogy, 
one  would  expect  convection  of  the  larger  particles  though  the  formation  of 
“fingers”  containing  high  concentrations  of  large  particles  extending  upward 
from  the  bottom  shear  plane.  If  a  similar  analysis  is  applied  to  the  viscous 
hydrodynamic  regime,  it  ma_>  provide  an  explanation  of  the  fingering  effects 
reported  by  Weiland,  et  al.16  Here,  particles  of  unequal  weight  and/or  size 
segregate  into  filaments  during  sedimentation.  Fluctuations  in  the  particle 
size  distribution  could  lead  to  more  rapid  sedimentation,  larger  shear  and  a 
horizontal  hydrodynamic  force  leading  to  further  separation. 

The  above  discussion  based  on  a  single  impurity  can  easily  be  extended 
to  describe  a  finite  concentration  of  impurities.  We  argued  that  an  approx¬ 
imate  solution  to  the  equation  of  motion,  Equation  ( B- 18),  for  shear  flow  in 
a  vertical  gravitational  field  is  p  ~  constant,  and,  hence.  /’  ~  T„  —  Gz.  where 
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T0  is  an  integration  constant.  In  equilibrium,  Equation  (3-21)  predicts  that 

DejjV2c  =  -VDejj  ■  VC, 

where 

DKff  oc  T*  &  (T0  -  gz)T 

Hence,  the  model  predicts  an  equilibrium  profile  of  the  form 

|  =  (4--0 

Several  considerations  should  be  taken  into  account  before  applying  these 
qualitative  conclusions  to  real  experiments.  First,  in  computing  the  diffusiv- 
ity,  we  have  used  the  Einstein  formula  which  assumes  viscous  flow  about  the 
large  particle  impurities.  However,  we  will  argue  in  the  next  section  that  the 
granular  flow  about  an  obstacle  is  likely  to  be  supersonic  and  contain  shock 
waves.  Thus,  it  is  not  clear  whether  the  effective  diffusivity  decreases  with 
increasing  R0  or  not,  as  the  Einstein  relation  would  predict.  We  raise  this 
qualification  because  the  microscopic  theory  would  predict  that,  impurities 
with  larger  radius  rise  more  rapidly,  whereas  the  argument  above  predicts 
the  opposite  (see  Equation  (3-20)).  One  possibility  is  that  both  predictions 
are  correct:  that  diffusivity  increases  with  particle  radius  if  the  impurity 
diameter  is  modestly  greater  than  the  grain  size,  but  the  opposite  holds  in 
the  limit  when  the  grain  s.ze  is  much  smaller  than  the  impurities.  In  this 
case,  it  should  be  very  interesting  to  study  experimentally  the  crossover  from 
microscopic  to  continuum-granular  behavior  as  a  function  of  impurity  size. 

Second,  we  note  that  experiments  which  are  reported  to  produce  the  seg¬ 
regation  effect  are  vert  ical  shaking  of  the  granular  material  rather  than  shear. 
It  may  be  that,  this  is  just  a  means  of  producing  an  increased  temperature 
and  higher  effective  diffusivity.  A  revealing  test  would  be  to  study  vertical 
shaking  versus  shear  and  to  discover  how  diffusivity  scales  with  particle  size. 
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Third,  the  above  estimate  will  break  down  if  the  temperature  change  due 
to  gravity  becomes  comparable  to  the  shear  induced  temperature: 


where  the  shear  rate  is  vo/ L.  If  we  use  a  grain  size  of  100  /im  in  a  box  of 
size  1  cm,  our  estimates  become  invalid  for  velocities  less  than  30  m/s ec  or 
so.  This  emphasizes  the  fact  that  the  local  shear  must  be  rather  large  for  it 
to  dominate  gravity  in  the  particle  dynamics. 

We  now  turn  to  the  issue  of  friction.  As  already  mentioned,  we  cannot 
have  infinitely  deep  granular  flow  regimes  because  the  temperature  rapidly 
decreases  away  from  the  boundary.  Even  in  the  case  where  the  temperature 
was  constant  (A  =  0),  including  gravity  will  cause  the  density  to  increase 
with  depth  and  limit  the  extent  of  the  granular  flow  layer. 

The  simplest  estimate  one  can  make  for  this  effect  is  to  take  an  all  or 
nothing  approach  to  frictional  forces.  That  is,  one  can  drop  friction  com¬ 
pletely  in  the  granular  flow  regime  but  require  that  the  flow  be  sufficiently 
vigorous  to  overcome  sliding.  This  means  that  the  ratio  of  shear  stress  r  to 
normal  stress  must  be  greater  than  tan0r,  where  cf)T  is  the  angle  of  repose. 
This  is  the  angle  that  a  sandpile  adopts  as  it  flows  (see  later).  Using  a  simple 
assumption  of  constant  density,  we  find  (Hanes  and  Inman1') 


T  =  T0 

p  =  Po-  pgfc  =  Po  +  pg(H  -  z) 
Ju 


so. 


H  -z  = 


IP9- 


This  predicts  that  the  layer  depth  for  the  system  should  increase  with  shear 
rate.  s;nce  r„  and  pG  increase.  If  we  measure  r0,  p„  and  the  layer  depth, 
however,  we  can  use  the  above  formula  to  find  an  effective  value  of  <t>r  which 
should  be  independent  of  shear  rate.  Experimentally,  this  appears  to  be 
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roughly  true,  but  there  is  a  systematic  increase  with  shear  rate  which  is  not 
explained  by  this  approach. 

How  might  one  do  better.  First,  we  have  not  properly  taken  into  account 
the  fact  that  the  density  is  not  constant.  Solving  for  the  density  involves 
finding  the  temperature  profile  (using  the  energy  equation  and  boundary 
conditions)  and  using  the  equation  of  state.  We  expect,  of  course,  for  the 
density  to  increase  as  we  go  downward  and  this  complicates  the  simple  linear 
relationship  given  above.  Since  this  effect  undoubtedly  depends  on  shear 
rate,  this  could  explain  the  deviation. 

Another  possibility  is  that  one  must  take  into  account  friction  for  the 
granular  flow  itself.  Roughly,  as  the  density  increases,  particle  contact  begin 
to  be  dominated  by  sliding  and  the  pressure  goes  up.  A  phenomological 
approach  to  this  idea  has  been  put  forth  by  Johnson  and  Jackson^13!  based 
on  critical  state  soil  theory,  but  this  approach  is  not  very  compelling.  Again, 
this  will  change  the  simple  relationship  between  depth  and  the  angle  of  repose 
and  could  explain  the  data. 

It  is  worth  pointing  out  that  one  must  be  careful  in  using  the  idea  of  a  fixed 
angle  of  repose.  As  has  been  stressed  by  Nageh18^  (and  is  apparently  well 
known  in  the  sand  community)  there  are  actually  two  angles  corresponding 
to  static  friction  and  dynamic  friction.  The  larger  angle  4>m  is  the  largest 
angle  for  a  sandpile;  once  this  angle  is  exceeded,  the  pile  flows  and  relaxes 
back  to  the  angle  of  repose  (f>T.  In  some  sense,  there  is  a  range  of  angles 
4>r  <  4>  <  <t>m  at  which  a  granular  fluid  can  coexist  with  a  granular  solid. 
The  natural  assumption  that  the  actual  equilibrium  point  lies  at  <j>r  may  not 
be  accurate.  This  too  may  offer  an  explanation  for  the  observed  variation  in 
the  friction  angle. 

There  is  a  completely  different  approach  to  granular  flow  that  starts  with 
the  static  granular  solid  and  attempts  to  model  flow  via  plasticity  equations. 
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This  appears  to  be  the  standard  methodology  employed  for  the  design  of 
grain  silos,  for  example.  We  would  like  to  point  out  that  there  is  increasing 
evidence  that  this  approach  fails;  this  evidence  includes  instabilities  of  the 
basic  assumed  profiles  (Schaeffer19)  and  the  apparent  lack  of  correlation  with 
experimental  findings  (see  Reference  2).  So,  although  there  is  as  yet  no  com¬ 
pletely  consistent  method  of  including  friction  in  the  granular  flow  theories, 
we  view  this  way  of  approaching  the  problem  as  more  productive  than  the 
idea  of  starting  with  the  solid  and  including  motion. 
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5  ACOUSTIC  MODES 


Let  us  return  to  the  simple  constant  temperature  profile  discussed  in 
Section  3.  As  Haff10  has  already  noted,  it  is  likely  that  there  will  be  important 
effects  due  to  density  variations  during  granular  flow.  This  is  true  because  the 
collision  frequency  and  hence  the  transport  depends  strongly  on  the  density. 
We  will  study  this  issue  by  focussing  on  the  acoustic  modes  of  the  granular 
fluid. 

Let  us  consider  the  linearized  version  of  the  fluid  equation;  we  will  assume 
that  the  acoustic  mode  is  propagating  in  the  y  direction,  normal  to  the  x-z 
shearing  plane.  The  equations  are 


(JL  +  L1L)  _  r<»>  (jL  +  UL) 

(pl°)  2  T(°)J  Vp<°> 


Here  the  superscript  means  the  value  in  the  base  solution  including  the 
density  and  temperature  dependence.  For  example 

p(°)  =  p(o)j<(o)Zl2 


In  the  last  equation,  we  can  simplify  the  right  hand  side  by  noting  the 
A  =  0  equilibrium  condition 


_  p(°) 


This  then  implies  that  the  last  two  terms  collapse  to  the  simpler  form 

pHyv 
Tl°)  ' 
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To  proceed,  we  take  a  time  derivative  of  the  density  equation  and  use  the 
momentum  equation  for  the  y  velocity  it'.  This  yields 


gv  c^y  .  .  .  ay 

dt2  ‘dy2  "  dy2  p1'^  '  1  dy2dt 

where  the  isothermal  sound  speed  c2  =  ^|r  and  c2  =  J^|  .  The  energy 
equation  using  the  above  simplification  has  the  form 

3  Mdr  Ao)d2r  p^dp'  r^r 

2P  dt  ~K  dy2  +  T(°)  ■ 


Assume  the  dependence  for  p\  T'  of  e^tky  ut\  we  find  the  eigenvalue  equation 

c2k 2  +  uj2  —  iuk2v'  c2k 2 


0  =  det 


2  P(o)X 


=**  +  k2  +  ^ 

X  X 


where 


X 


u<°> 

A :<•> 

3/2p(°) 

p(°) 


This  gives  us  a  fourth  order  equation  for  uj(k)  of  the  form 


1  - 


lUJV 


+  k2 


—UJ 


UJr 


+  —  1  - 


—lUJV 


—  IUJ 


1 


IUJV  \  —luJC 


UJ 


3  (o)  "b  2  ^ c } 

jr  °  A  \Ct 


The  remainder  of  this  section  will  be  devoted  to  examining  the  consequences 
of  the  above  dispersion  relation. 


Before  proceeding  it  is  worthwhile  to  digress  and  consider  the  same  prob¬ 
lem  for  normal  fluid.  There,  ujc  —  0  since  there  is  no  dissipation.  For  gases 
the  kinematic  viscosity  v'  is  negligible  compared  to  the  thermal  diffusivity. 
If  we  define  a  frequency  u;x  =  c2/\  we  can  distinguish  two  limits 
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(a)  u  <u;x: 


c^2  (1  +  fife)  =  u;2 

=4*a 


where  c3  is  the  adiabatic  sound  speed.  There  is  also  an  absorption  coefficient 
where 


u> 


k  = - b  iau2,  a  ~  x 

c. 


(b)  u>  >  cjx! 

Now  we  get  approximately  isothermal  propagation  since  the  temperature 

field  is  frozen  on  the  time  side  of  the  density  fluctuation.  Now 

,  u>  ...  1 

k  — - b  la  a  ~  . 

ct  X 

Of  course,  u'  is  not  exactly  zero  and  eventually  becomes  important.  If  u>  is 
much  bigger  than  w 'u  =  c2Ju',  we  then  get  a  balance  of  the  form 

—iu>kA 

U)u' 

For  air,  this  occurs  at  frequencies  of  109  Hz;  for  water  at  1012  Hz.  Hence  these 
viscous  regimes  are  not  of  much  practical  importance.  There  is  of  course  also 
a  second  mode  at  small  frequency 


which  corresponds  to  thermal  diffusion,  with  no  velocity  or  density  pertur¬ 
bation. 

We  now  return  to  the  granular  fluid.  The  major  difference  is  that  tem¬ 
perature  fluctuations  are  damped  even  at  low  frequency.  Assuming  that  u  is 
smaller  than  all  other  frequency  (ujc,ljx,  and  uv  =  c2/i/')  we  obtain 

,  U?  2  ,  / 

k  — - bzua  a  ~  1  /a )c. 

c, 
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That  is,  the  propagation  is  isothermal!  At  high  frequencies,  we  again  have 
purely  damped  motion  with  two  possible  modes 


Again,  there  is  also  a  second  mode  at  small  frequency  which  is  purely  damped, 
as  expected  from  the  discussion  above. 

Why  are  these  modes  interesting?  First  of  all,  it  offers  a  completely  new 
methodology  for  testing  the  continuum  theory  of  granular  flow.  It  is  certainly 
possible  to  set  up  a  shear  flow  and  measure  all  the  relevant  material  properties 
(It  is  also  possible  to  calculate  them  within  the  context  of  a  kinetic  theory 
which  starts  from  particle  collisions,  but  this  is  not  likely  to  yield  much 
quantitatively  valid  information).  Once  all  these  parameters  are  known,  we 
would  have  a  prediction  of  the  sound  speed  and  the  relationship  between 
density  and  temperature  fluctuations.  Via  exciting  these  modes  with,  say,  an 
oscillatory  variation  of  the  driving  plate,  one  could  test  this  entire  framework. 

There  is  a  second  reason  why  sound  modes  are  vital.  To  see  this,  recall 
that  in  a  granular  flow  all  velocities  are  ultimately  related  to  the  shear  rate 
imposed  on  the  flow.  For  example,  in  Section  3  we  saw  that  the  typical 
fluctuation  velocity  vj  ~  T1^2  is  proportional  to  7 a  where  a  is  the  particle 
size.  Let  us  therefore  consider  the  Mach  number  of  a  typical  flow 

M  —  vfct. 

Now  v  ~  7 L  where  L  is  the  system  depth  and  ct  ~  ~  T1^2  ~  7(7.  Hence 

the  ratio  A/  ~  L/cr  times  some  complicated  function  of  density.  Since  L 
must  be  at  least  tens  of  particle  sizes  for  a  macroscopic  theory  to  make  any 
sense,  we  have  the  simple  result  that  many  typical  flows  will  be  supersonic. 

The  above  statement  again  has  many  interesting  consequences  for  testing 
the  validity  of  continuum  theories.  We  would  predict  for  example  that  placing 
an  obstacle  in  the  flow  would  cause  the  formation  of  shock- like  fronts  in  the 
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density.  It  should  be  noted,  however,  that  the  highly  dissipative  nature  of 
granular  fluids  (i.e.  large  viscosity  and  temperature  fluctuation  dumping)  will 
tend  to  significantly  broaden  such  fronts.  Also,  any  nascent  instability  in  the 
steady  flow  would  couple  to  compressibility  in  a  highly  nontrivial  manner. 
Such  effects  form  the  subject  matter  of  the  next  section. 

Although  most  of  the  results  in  this  work  have  dealt  with  the  granular 
flow  system,  we  would  like  to  remark  on  density  fluctuations  in  a  viscous 
suspension.  As  mentioned  at  the  end  of  Section  2,  one  proposed  model  for 
viscous  suspensions  would  be  a  similar  set  of  hydrodynamic  equations  with 
the  replacement 

p  oJT 

Doing  this  makes  the  pressure  independent  of  the  particle  density.  This 
removes  the  basic  driving  force  for  wave  phenomena  and  one  can  show  that 
all  modes  are  damped.  The  two  modes  respectively  are  temperature  which 
is  just  damped  and  density  which  is  diffusive,  i.e.  u  ~  ik2.  This  is  quite 
reasonable  since  the  only  sound  waves  expected  in  the  system  rely  on  the 
compressibility  of  the  ambient  fluid  to  transmit  fluctuations.  This  is  not 
irrelevant  in  lubrication  force  dominated  regimes. 


6  INSTABILITIES  IN  GRANULAR  FLOW 


In  this  section,  we  consider  possible  instabilities  which  may  transform 
steady  granular  flow  patterns  into  more  complex  spatial  and  temporal  struc¬ 
tures.  Again,  a  major  motivation  is  to  provide  a  variety  of  qualitative  tests  of 
the  idea  of  treating  solid-liquid  systems  via  a  one  phase  continuum  theory.  In 
fact,  there  is  some  evidence  of  anomalous  fluctuations  associated  with  “flow 
disorder”  in  experiments  of  Savage  and  Sayed  as  well  as  the  hopper  flow  of 
Baxter  et  al.  In  the  hopper  flow,  these  fluctuations  are  probably  connected 
with  observed  density  waves  which  propagate  through  the  system. 


Let  us  focus  on  a  Taylor-Coutte  system  as  an  example  of  a  granular  flow 
instability.  First,  let  us  recall  the  classical  case  of  an  incompressible  flow. 
Our  base  flow  is  purely  azimuthal 

V(o)  =  u(r)0. 


Assuming  a  small  perturbation  vr,  i>@,p,  we  have 

- 1 


p^dvz 

dt 

P(o)dv  e 
dt 


dp 

"d~z 

-Vr 


d 


r  V0)^(r3n(r)). 


Now  the  vorticity  of  the  original  profile  Q(r)  =  u©(r)/r  is  arbitrary  at  zero 
viscosity.  If  we  include  viscosity  in  the  base  solution  while  still  neglecting  it 
for  the  perturbative  analysis,  we  should  take 


n  =  0,  +  % 

rl 


since  the  original  momentum  equation  takes  the  form 

/  d2  1  d  1 


Returning  to  the  stability  calculation,  incompressibility  demands 


kzVz  =  -krvr. 

We  can  eliminate  by  recognizing  that 

1  dp  1  kT  dp  —kT  dvz 


kz 


dvr 
~dt  ' 


p(°)  dr  p(°)  kz  dz  kz  dt 
Substituting  this  and  also  the  equation  v&  into  the  equation  for  dvr/dt,  we 
get  the  dispersion  relation 


2th  d  .  .  . 

ilo  =  - 7 — (r  s l)  +  iu> 

u>  dr 


V 

k. 


or 


LO  = 


k2z<t> 
k2  +  k] 


^  2nd  t  2 
<t>  =  — "H). 

r  dr 

The  instability  criteria  thus  reduces  to  the  well-known  Rayleigh  criterion 
<f>  <  0  for  all  r.  This  can  be  employed  to  show  that  as  Taylor-Couette  cell 
with  counter-rotating  cylinders  is  always  unstable  whereas  stability  persists 
in  the  co-rotating  case  until  the  rotation  ratio  exceeds  a  critical  value. 


We  now  would  like  to  understand  what  happens  if  one  reconsiders  the 
above  analysis  in  the  case  of  a  granular  fluid.  The  first  point  is  that  there 
may  be  significant  stratification  due  to  the  centrifugal  force.  Specifically, 
momentum  balance  in  the  r  direction  requires 

dp  _  pv2 
dr  r 

and  using  p  ~  pT 

rr  r  [rr'n2(r’) 

PT  =  P^cW  J-^dr 

where  V'(r)  =  rfl(r).  The  corresponding  vorticity  profile  now  obeys  the 
modified  equation 

n  (  1)2  9p  dn 
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whose  solution,  using  p  =  p0pTll2,  is 

ft  =  ft0  +  fti  f  -4 pT1/2{r ')  dr'. 

J  T  o  7* 

The  system  of  equations  for  pT  and  ft  is  now  closed  by  writing  down  the 
energy  equation.  Clearly,  this  is  quite  complicated.  However,  our  main  point 
is  that  there  may  be  a  strongly  inhomogeneous  distribution  of  density  and 
temperature  in  anything  other  than  the  thin  gap  limit. 


for  density  variations,  and  the  energy  equation,  in  the  limit  u>  »  wc  (and  all 
other  dissipative  frequencies)  has  the  form 


df  .  dTW  2p(°)~  A 

dt  +Vr  dr  ~  37^)V'U' 

The  key  difference  between  these  equations  and  the  ones  from  the  normal 
fluid  is  the  presence  of  density  fluctuations  which  then  couple  to  temperature 
fluctuations. 


We  will,  for  simplicity,  deal  with  the  highly  supersonic  limit 

u>  >  kct. 
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In  this  limit,  we  can  drop  the  V  •  V  term  in  the  temperature  and  density 
equation,  since  convection  is  much  more  rapid  than  transmission  by  acoustic 
waves.  So 


Vo  = 
V, 


VT  6  2 


iur  dr 

2QVo  + 


=  ±( 
ILJ  \ 


m_\ 

Vp(°) ) 


+  0 


'ctk' 

U  , 


From  the  density  equation 


P  = 


Vr  dpV) 


iu  dr 

substituting  the  Vo  equation  to  the  VT  one,  we  find  the  dispersion  relation 

2  ve  dpV)  n(ctk\ 

“  _$  +  ^w;“8r+0(^r)' 


Comparing  this  formula  to  the  preceding  one,  we  see  the  presence  of  a 
nev>  term  related  to  the  stratification  If  the  material  is  more  dense 

at  larger  radius  (certainly  the  expected  effect  of  the  centrifugal  term),  this 
will  tend  to  stabilize  the  system,  essentially  via  a  Rayleigh-Taylor  density 
interchange.  The  upshot  is  that  we  observe  an  interesting  coupling  between 
angular  momentum  and  density  stratification,  which  will  change  the  stability 
criterion  as  well  as  the  character  of  the  unstable  mode. 

The  above  analysis  suggests  that  it  would  be  very  interesting  to  systemat¬ 
ically  investigate  all  the  “typical”  fluid  instabilities  for  granular  liquids.  For 
the  completely  supersonic  resinc  lo  kct  (as  well  all  dissipative  frequencies), 
this  amounts  to  including  compressibility.  For  smaller  u>  (and  hence  for  the 
onset  of  all  non-oscillatory  instabilities)  the  strong  damping  of  temperature 
fluctuations  will  be  important  (and  presumably  ensure  that  growing  modes 
are  mostly  isothermal),  as  will  the  coupling  to  density  fluctuations.  Again, 
these  hydrodynamic  tests  of  the  granular  flow  concept  are  much  more  in¬ 
teresting  than  mere  viscometric  tests  which  have  already  been  performed  to 
date. 
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7  CONCLUSIONS 


In  this  report,  we  have  explored  the  single- phase  continuum  model  pro¬ 
posed  by  Jenkins  and  Savage  (1982)  to  describe  granular  flow.  Many  of  our 
results  can  be  extended  to  make  predictions  for  dense  suspensions,  which 
can  be  described  using  a  similar  theoretical  framework.  The  goal  has  been  to 
develop  clear,  qualitative  predictions  that  appear  to  be  amenable  to  exper¬ 
imental  verification.  In  particular,  our  analysis  suggests  that  an  increased 
battery  of  hydrodynamic  tests  on  granular  (and  dense  suspension)  systems 
could  be  developed  that  would  be  more  effective  in  sorting  out  valid  model  ap¬ 
proaches  compared  with  proceeding  directly  to  expensive,  quantitative  tests. 
We  would  strongly  oppose  large-scale  simulations  or  quantitative  tests  of  this 
or  any  other  model  until  the  model  passes  these  types  of  qualitative  tests  and 
its  assumptions  have  been  validated. 

The  key  results  from  our  analysis  of  granular  flow  based  on  the  single¬ 
phase  continuum  model  include  the  following. 

•  We  have  predicted  the  flow  profiles  for  simple  boundary  conditions, 
such  as  plane  Couette  flow  and  channel  flow.  From  these  studies, 
macroscopic  and  microscopic  criteria  for  the  model’s  validity  have  been 
developed.  For  the  potentially  valid  range  of  concentration  and  veloc¬ 
ity,  pressure  and  shear  stress  should  scale  quadratically  with  shear  rate, 
and  there  should  be  a  near-balance  between  viscous  heating  and  colli- 
sional  dissipation  on  microscopic  scales. 

•  We  have  argued  that  the  model  predicts  segregation  effects  when  there 
are  two  (or  more)  particle  sizes  considered.  In  shear  flow,  large  particles 
should  rise  in  a  gravitational  field  compared  to  small  grains,  perhaps 
producing  double-diffusive  fingering  in  some  limits. 
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•  We  have  derived  the  acoust'c  dispersion  relations  and  shown  that  they 
predict  granular  flow  to  be  supersonic.  This  result  can  be  verified  by 
measuring  the  sound  speed  in  shear  flows.  Also,  the  flow  past  bodies  in 
granular  shear  flows  should  reveal  disturbances  only  past  the  leading 
edge  of  the  body  and  broadened  shock  fronts  characteristic  of  dissipa¬ 
tive  flow. 

•  We  point  out  that  the  analysis  of  flow  instabilities  is  complicated  by  the 
finite  compressibility  of  the  solid-fluid  mixture.  For  example,  the  finite 
compressibility  leads  to  the  prediction  of  density  interchange  (Rayleigh- 
Taylor  instabilities)  in  addition  the  usual  angular  momentum  inter¬ 
change  in  standard  (cylindrical)  Couette  flow. 


We  believe  that  these  predictions  suggest  a  series  of  experiments  that  can 
form  the  basis  of  an  exciting,  insightful,  and  (hopefully)  inexpensive  program 
of  research  on  solid-fluid  flow.  Because  of  their  qualitative  nature,  these  tests 
should  provide  an  empirical  verification  of  theoretical  models  that  is  less 
sensitive  to  the  details  of  the  system.  In  the  Appendix,  we  outline  existing 
and  potential  approaches  for  imaging  multiphase  flow,  which  would  be  a  key 
component  in  developing  the  experimental  program. 
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A  IMAGING  SOLID-GAS  AND  SOLID-LIQUID 
FLOWS 


In  Section  I,  we  commented  on  the  need  for  an  expanded  range  of  qualita¬ 
tive  hydrodynamic  tests  of  theories  before  considering  expensive,  large-scale 
computer  simulations.  For  example,  several  different  types  of  experiments 
were  suggested  in  the  discussion  of  the  single- phase  continuum  model  of 
solid-fluid  flow.  A  common  need  in  all  of  these  qualitative  tests  is  detailed 
information  about  particle-fluid  positions  and  velocities. 

Imaging  techniques  of  various  types  have  been  developed  by  many  re¬ 
searchers  to  gather  this  information.  Some  of  these  approaches  are  restricted 
to  special,  ideal  test  systems;  others  can  be  applied  in  the  real  plant  envi¬ 
ronment,  as  well.  However,  as  we  discuss  below,  none  of  these  techniques  is 
sufficiently  powerful  to  provide  a  general  solution  to  the  imaging  problem. 
Thus,  an  important  component  of  experimental  research  in  this  area  should 
be  the  improvement  of  imaging  techniques. 

Because  this  issue  is  so  closely  tied  to  the  development  of  theory,  we 
thought  it  might  be  useful  to  comment  on  present  and  projected  imaging 
methods.  Below,  we  briefly  discuss  imaging  techniques  for  both  solid-liquid 
and  solid-gas  flows.  The  flows  important  for  applications  often  have  a  rela¬ 
tively  rapid  time  development.  An  ideal  imaging  system  would  provide  two- 
or  three-dimensional  images  of  the  different  phases  in  the  flow  with  good 
time  resolution  and  with  little  or  no  perturbation  of  the  flow.  One  must  also 
bear  in  mind  that  the  flow  often  creates  a  hostile  environment  for  sensors, 
for  example,  by  sand-blasting.  As  we  will  find  below,  these  goals  are  difficult 
to  meet. 
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A.l  Non-Invasive  Techniques 


Electrical 


Capacitance  probe:  In  this  technique,  used  by  Louge  and  collaborators 
at  Cornell,  a  capacitive  probe  is  located  along  the  wall  of  a  pipe  contain¬ 
ing  a  particulate  suspension.  The  probe  is  sensitive  to  changes  in  dielectric 
constant  between  the  fluid  and  solid.  With  care,  excellent  sensitivity  can 
be  achieved  using  commercially  available  electronics.  This  technique  can 
be  readily  extended  to  arrays  of  capacitive  probes  to  provide  some  sp  > ,  J 
information. 

Present  geometry  by  Louge  et  al.  needs  to  be  improved  for  better  control 
of  fields.  ADVANTAGES:  fast;  sensitive;  rugged  probe;  relative  cheap.  DIS¬ 
ADVANTAGES:  local  probe  only  senses  flow  near  surface;  pickup  area  not 
sharply  defined. 


Inductive  pickup:  A  pickup  coil  is  wound  around  the  outside  of  the  pipe. 
For  particular  flows,  the  particles  often  acquire  an  electrostatic  charge,  which 
produces  a  pulse  in  the  pickup  coil  when  the  particle  pass  through.  One 
can  also  imagine  charging  particles  with  an  electrical  discharge.  Multiple 
pickups  could  be  used  along  the  pipe.  ADVANTAGES:  fast,  cheap;  rugged. 
DISADVANTAGES:  shielded  by  metal  pipes;  pickup  has  spatial  resolution 
only  along  the  pipe;  works  only  for  charged  particle  such  as  sand  grains  in 
gases;  details  of  charging  difficult  to  control  and  may  influence  flow. 

Acoustic 


Sound-passive  microphones:  Conventional  contact  microphones  are  at¬ 
tached  along  the  outside  of  the  pipe.  Spectrum  analysis  may  provide  some 
useful  information  about  the  flow.  ADVANTAGES:  fast;  cheap;  portable; 
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could  be  installed  along  pipelines  in  plants  as  an  online  diagnostic.  DISAD¬ 
VANTAGES:  nonlocal  pickup;  hard  to  analyze  quantitatively. 

Pressure  probes:  Conventional  pickups  are  used  to  measure  the  pres¬ 
sure  in  the  flow.  For  fluidized  beds,  this  is  a  good  measure  of  the  average 
particulate  content.  ADVANTAGES:  reliable;  cheap;  simple  interpretation. 
DISADVANTAGES:  poor  spatial  resolution;  only  work  near  pipe  surface. 

Optical 


Optical  imaging  with  tracer  particles:  In  this  technique,  visible  tracer 
particles  are  added  to  a  multiphase  flow  which  is  made  transparent  by  us¬ 
ing  transparent  liquid  and  particles  with  matched  indices  of  refraction.  If 
the  tracers  have  properties  similar  to  particles,  then  one  images  the  particle 
flow;  if  they  are  flow  visualizers  such  as  “fish  scales,”  then  one  images  the 
fluid.  These  techniques  and  variants  are  used  for  the  study  of  ideal  systems 
in  the  laboratory  (e.g.,  J.  R.  Kadambi,  Case  Western  Reserve).  ADVAN¬ 
TAGES:  visual;  3D,  fast;  easy  to  interpret;  uses  conventional  photographic 
and  computer  techniques;  relatively  low  cost.  DISADVANTAGES:  limited 
to  ideal  systems  only. 

Optical  imaging  with  laser  induced  fluorescence  or  color  change:  This 
technique  also  requires  transparent  multiphase  flow,  and  hence  is  limited  to 
ideal  systems.  As  in  the  work  of  R.  Falco  (Michigan  State),  a  high  power 
laser  is  used  to  induce  a  short-lived  fluorescence  or  color  change  in  either  the 
particles  or  the  liquid  of  a  multiphase  flow  in  a  well-defined  pattern,  such  as 
a  grid.  Thus,  the  flow  is  tagged  at  one  instant  with  a  pattern.  The  distortion 
of  this  pattern  in  time  allows  one  to  derive  the  local  motion  of  the  particles 
or  fluid,  including  the  local  vorticity.  We  strongly  recommend  that  the  Falco 
and  Kadambi  (previous  technique)  groups  be  in  close  contact  to  share  tech¬ 
nology  and  speed  development  of  the  optimal  technique.  ADVANTAGES: 
visual;  3D;  fast;  easy  to  interpret;  uses  conventional  photographic  techniques; 
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computer  image  processing  permits  mapping  of  local  velocity  and  vorticity 
profiles.  DISADVANTAGES:  precise  index  match  limits  technique  to  ideal 
systems;  relatively  high  cost  of  laser  and  computer  image  processor. 

Medical  Techniques 

In  developing  the  ideal  imaging  technique  for  solid-fluid  flow,  one  is  faced 
with  a  problem  already  encountered  in  medicine,  i.e.,  how  to  image  solids 
suspended  in  water  using  non-invasive  techniques.  An  enormous  effort  has 
been  made  to  develop  these  techniques  in  the  medical  area,  with  many  no¬ 
table  successes:  sophisticated  x-ray  imaging,  x-ray  tomography,  ultrasonic 
and  magnetic  resonance  imaging.  The  equipment  required  is  often  very  ex¬ 
pensive  (over  1M  machine),  but  is  readily  available  in  many  medical  centers 
located  throughout  the  country  for  use  in  cooperative  research  or  simply  for 
an  hourly  fee.  For  example,  Behringer  cite  BEH  has  collaborated  with  the 
Medical  School  at  Duke  University  to  make  time-resolved  x-ray  images  of 
flowing  sand.  The  use  of  these  facilities  can  be  a  sensible  and  cost-effective 
solution  for  those  experiments  which  can  be  transported. 

X-ray  imaging:  Standard  x-ray  cameras  can  be  used  to  image  multiphase 
flows.  Digital  image  subtraction  can  be  used  to  detect  changes  in  flow,  as 
in  the  research  of  Behringer  (Duke).  ADVANTAGES:  standard  technique 
easily  accessible  in  hospitals;  good  time  resolution  for  short  bursts  of  images. 
DISADVANTAGES:  2D  image;  expensive  if  not  in  hospital;  low  contrast  if 
used  for  liquid-gas  mixtures. 

X-ray  tomography.  Conventional  CAT  scanner.  ADVANTAGES:  3D  im¬ 
age;  available  at  many  medical  centers;  excellent  detail.  DISADVANTAGES: 
too  slow  for  many  applications;  prohibitively  expensive  outside  of  hospital. 

NMR  imaging:  Conventional  magnetic  resonance  imager.  ADVANTAGES: 
3D  image;  available  at  many  medical  centers;  excellent  detail  and  contrast. 
DISADVANTAGES:  very  slow;  prohibitively  expensive  outside  of  hospital. 
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Ultrasonic  imaging :  Conventional  medical  ultrasonic  imaging  instrumen¬ 
tation.  The  application  to  multiphase  flows  is  questionable  but  worth  inves¬ 
tigation.  Its  success  will  depend  upon  the  acoustic  background  created  by 
the  flow  at  ultrasonic  frequencies  and  by  clutter  created  by  multiple  scattered 
returns.  ADVANTAGES:  available  at  hospitals;  portable  apparatus;  mod¬ 
erate  expense.  DISADVANTAGES:  less  resolution;  slow;  high  background 
noise  level  and  clutter  from  flow;  requires  liquid  medium. 

A. 2  Invasive  Techniques 

Capacitance  probes:  The  technique  is  similar  to  the  non-invasive  counter¬ 
part  above,  except  that  the  probe  is  positioned  in  the  flow.  Arrays  of  probes 
can  also  be  used  with  greater  perturbation  of  the  flow.  If  the  prediction  that 
the  flow  is  supersonic  is  correct  (Section  4),  then  probe  only  disturbs  flow 
past  the  measuring  point.  ADVANTAGES:  as  above;  mapping  of  concentra¬ 
tion  profiles  possible;  probes  can  be  small  and  rugged.  DISADVANTAGES: 
single  point  measurements;  pickup  area  not  sharply  defined;  perturbs  flow. 

Velocimeter  probes:  Laser  induced  fluorescence  can  be  used  to  tag  a  single 
particle  at  a  given  instant.  The  motion  of  the  particle  is  then  tracked  using 
optical  probes  (for  example,  an  array  of  fiber-optic  pickups  or  photodiodes), 
and  the  timing  of  its  passage  used  to  determine  the  magnitude  and  direction 
of  its  velocity.  The  laser  light  can  be  brought  to  the  end  of  a  small  probe  using 
fiber-optics.  ADVANTAGES:  mapping  of  velocity  profiles;  fast;  potentially 
small  and  rugged.  DISADVANTAGES:  single  point  measurement;  perturbs 
flow. 

Sampling  probes:  Here  a  hollow  probe  is  simply  used  to  collect  a  sample 
of  the  particles  at  a  given  position.  ADVANTAGES:  mapping  of  particle 
characteristics  vs.  position  to  test  for  segregation,  etc.  DISADVANTAGES: 
slow;  time  averaged;  perturbs  flow. 
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